Superradiance and phase multistability in circuit quantum electrodynamics 
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By modeling the coupling of multiple superconducting qubits to a single cavity in the circuit- 
quantum electrodynamics (QED) framework we find that it should be possible to observe super- 
radiance and phase multistability using currently available technology. Due to the exceptionally 
large couplings present in circuit-QED we predict that superradiant microwave pulses should be 
observable with only a very small number of qubits (just three or four), in the presence of energy 
relaxation and non-uniform qubit-field coupling strengths. This paves the way for circuit-QED 
implementations of superradiant state readout and decoherence free subspace state encoding in sub- 
radiant states. The system considered here also exhibits phase multistability when driven with large 
field amplitudes, and this effect may have applications for collective qubit readout and for quantum 
feedback protocols. 
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I. INTRODUCTION 

In 1954, Dicke [1 noted that the intensity of radia- 
tion from a gas could be enhanced by confining it in 
a region smaller than the wavelength of emitted radia- 
tion, a phenomenon called superradiance. The collective 
emission intensity, which scales with ex N 2 , is of funda- 
mental interest in quantum optics and has been the sub- 
ject of many theoretical and experimental works [2HI]. 
With greater control over a range of quantum systems, 
it is now possible to observe superradiance in quantum 
dots [5], Bose Einstein condensates [5] and nitrogen va- 
cancy centers in diamond [7] . Furthermore superradiance 
and the related phenomena of subradiance have quantum 
information based applications such as decoherence free 
subspace state encoding in subradiant states [H] and su- 
perradiant state readout [9]. 

Superradiance has been difficult to observe experi- 
mentally as decoherence and dissipation destroy the re- 
quired build up of correlations within the ensemble [2H1] . 
Plagued by these losses, experiments have required large 
ensembles to observe the ex TV 2 intensity characteristic 
of superradiance [SHU HMT^] . This has inhibited any de- 
tailed study of emission dynamics or measurement of cor- 
relations within the ensemble, leaving the phenomenon 
demonstrated in principle yet not explored in detail. 

Another phenomenon that is difficult to observe ex- 
perimentally is phase bistability. In 1991, Alsing and 
Carmichael, and independently Kilin and Krinitskaya, 
showed that a strongly driven qubit-cavity system can 
undergo a symmetry breaking phase transition to a state 
where the intracavity field is displaced in phase, depen- 
dent on the qubit's state [131 Q3]. This state has the 
remarkable property that the qubit-cavity system can be 
described as two coupled, driven and damped harmonic 
oscillators. Kilin and Krinitskaya also predicted that a 
similar phenomenon, phase multistability, occurs when 
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FIG. 1: Five transmon qubits (dark blue) coupled to the 
quantized field of a transmission line resonator (sinusoids). 



many qubits couple to a single cavity mode [15] , Phase 
multistability describes the phenomenon, where a coher- 
ently driven multi-qubit cavity system displaces the cav- 
ity field in phase depending on the collective state of the 
qubits. Phase multistability for N > 1 has yet to be 
observed. 

Recently, several promising experiments in circuit- 
QED have demonstrated various quantum information 
based tasks such as three qubit entanglement, multi- 
qubit measurement and the realisation of simple quan- 
tum algorithms [TfjrtTB] . With such a large degree of con- 
trol over system dynamics the question naturally arises, 
"Can we use superconducting circuits to probe many 
body collective effects that were previously inaccessible 
to experiment?" . By demonstrating that small sample 
superradiance and phase multistability can be observed 
in circuit-QED, we provide an affirmative answer to this 
question. Circuit-QED is an ideal platform to observe 
small ensemble superradiance and phase multistability 
due to small qubit losses and large qubit-field coupling 
rates [19]. 

This article is organized as follows. In section [TTJ we 
show that circuit-QED can probe the small ensemble 
regime of superradiance, allowing its detailed investiga- 
tion. In section |III| we extend the system by adding a 



coherent driving of the resonator mode. This leads to 
the phenomenon of phase multistability which is anal- 
ysed in detail. Our conclusions are presented in section 

nvi 



II. SUPERRADIANCE 

Here, we study several transmon qubits coupled to the 
field of a transmission line resonator (TLR) (Fig. [l} . In 
the frame rotating at angular frequency u> the dynamics 
of the system are described by the master equation, 

N s p 

p = -i[H, P \ + * z v[a]p + Y, 1 i v ^p+ 1 i v ^m 

where 7* and 7^ are for the jth qubit the energy relax- 
ation and dephasing rate respectively, k is the resonator 
decay rate, T>[A]p = 2ApA^ - A* Ap - pA^A and h = 1. 
The qubit-resonator system evolves in the interaction pic- 
ture under the Tavis-Cummings Hamiltonian [20j, 
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where, the jth qubit couples to the field at the rate gN,j, 
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u) is the detuning from the jth qubit tran- 



sition frequency w gj - and A r = w r — u is the detuning 
from the resonator frequency u r . The system is initially 
taken to be in the state where the qubits are prepared 
in a symmetrically excited state and there are no pho- 
tons in the cavity, |^(0)) = |0) ® |e, e, . . . , e). Due to the 
Tavis-Cummings interaction, energy is transfered to the 
resonator. We probe the superradiant emission of these 
photons from the TLR in the "bad cavity" limit [21], 
where the photons escape the cavity before re-absorption 
by the qubits. 

In order to obtain an analytic expression for the inten- 
sity of emitted photons from the TLR, it is necessary to 
make several assumptions. First, we assume that each 
qubit couples identically to the field mode (gN,j ~ §n), 
as non-uniform coupling rates induce subradiant tran- 
sitions. This can be achieved by optimizing the place- 
ment of qubits at field antinodes [TH1 [TU] . Secondly, in 
the bad cavity limit it is assumed that the field in the 
resonator decays at a faster rate than the average qubit- 
field coupling rate, k ^$> cjn- Lastly, it is assumed that 
the dissipative processes, dephasing (7?) and energy re- 
laxation (7?), are negligible on the time scale of interest, 

k>3at >7j,7f 

As the resonator field decays on a time scale much 
faster than that of the qubits we can adiabatically elim- 
inate the field degrees of freedom [22 . The master 
equation for the reduced density operator describing the 
qubits only is given by, 
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where, T = n/2 + iA r , and it is assumed that the res- 
onator is initially empty. The collective operators are 
defined as, Ji — Y^j=i a ji where i = {+,— ,z}, and erj 
denote the individual qubit Pauli matrices. Equation 
(pj) corresponds to the superradiance master equation 
(SRME) that was derived for atomic ensembles, after a 
suitable parameter substitution [2] [23] . 

We now seek an expression for the intensity of photons 
escaping the TLR, 7jv(i). Expanding fcfy in the Dicke 
basis PQ, the probability that the system is in one of 
the Dicke states \l,m) is, P(l,m,t) = (l,m\p q (t)\l,m}, 
where, J 2 \l,m) = 1(1 + l)\l,m), J z \l,m) = 2m\l,m) and 
N/2 > I > \m\ > 0. Using the SRME ^ the population 
rates follow [24], 
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P(l,m,t) = r^r[(l-m)(l + m + l)P(l,m + l,t) 

\ T \ 

- (l + m)(l-m + l)P{l,m,t)]. 
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Here, we consider the initial condition that all qubits are 
excited, i.e., P(l,m,0) = P(N/2, N/2,0), although this 
approach can easily be extended to take into account 
other initial conditions [5S]. As I is conserved by the 
SRME we introduce the variable, n = I — rn = 0, 1, ..., N, 
which corresponds to the number of photons emitted 
from the resonator when the system is in the initial state 
m = N/2. Equation Q can now be rewritten, 

P(n,T) = (N-n + l)nP(n-l,T)-{N-n)(n + l)P{n, 

where we have rescaled the time t = jt and 7 = 
Kgff/\T\ 2 . To proceed, we Laplace transform (pM), subject 
to the full excitation initial condition, P(n,0) = 5 n ^, 

sQ(n,s) — S n fl = (N — n + l)nQ(n — 1, s) 

- (N-n)(n+l)Q{n,s), (6) 

where, Q(n, s) is the Laplace transform of P(n, r). From 
ffib we find Q(0,s) = l/(.s + N). Continuing recursively 
we find, 
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Upon inverting the transform ¥fh we obtain the popula- 
tions of the Dicke states, P(ti,t) (23] ■ The intensity of 
photons emitted from the TLR can be found from the 
Dicke state populations using, In( t ) = ■$- ^2 n n P( n , T )- 
The system is superradiant if the maximum intensity, 
I^ ax , is greater than N, the initial intensity of ./V inde- 
pendent qubits pp. 

Solving the SRME for the intensity profiles for exper- 
imentally relevant systems of TLRs with ./V = 3, 4 or 5 
qubits we find, 



I 3 (t) = 3e- 3T (12r-7) + 24e- 4T . 

I 4 (t) = (72r + 96)e- 6T +4e~ 4T (36r-23) 

I 5 ( T ) = |[162e- 9T + 16e- 8r (24r-l) + e- 
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The system is superradiant if the maximum intensity, 
r^ ax , is greater than TV, the initial intensity of N inde- 



pendent qubits 0. As I^ ax 
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FIG. 2: (Color Online) Intensity of radiation from the res- 
onator Iff(Rit) obtained from numerical solution of the com- 
plete master equation, where Ri = r(A r = 0) = Ag% / k,. 
(a), N - 3. (b), N = 4. (c), N - 5. Superradiant 
emission is colored red. Parameters used are gz^fl-n = 
(83.7,85.7,85.1) MHz [B], g 4 ,j/2w = (69.4,69.1,68.6,69.7) 
MHz, gs,j/2-K = (59.0,59.4,59.9,60.9,60.7) MHz, and 
(K,7l,7|)/27r= (2000,0.19,0) MHz. 



and, I™ ax ss 6.9 > 5, these systems exhibit superradi- 
ance. 

In Fig. [2] the intensity of emission from the ensembles 
is shown for a range of experimentally feasible parame- 
ters. It is clear that superradiance (colored red) can be 
observed over a large range of cavity decay rates, n. The 
approximate solutions pa|c) closely resemble the results 
in Fig. [2] for k >• g N ^> 7?, 7?- For each system the 
superradiant peak is large enough to be easily resolved 
using existing detection schemes [2l)H2"5] . 

The above results demonstrate that it is possible to 
observe small sample superradiance in a TLR in the 
presence of energy relaxation and non-uniform coupling 
strengths. Superradiance can also be observed for TV > 5 
qubits, provided k^> (jn ^> Ijilp an d gN.j ~ 9n, how- 
ever this may be more difficult to achieve experimentally. 
For a TLR with resonance frequency cj r , the length de- 
termines the number of qubits that can be effectively 
coupled at the maximal coupling strength as each qubit 
is typically placed at the field antinodes. When N is in- 
creased the coupling rate to each qubit is reduced as the 
coupling rate §n is inversely related to the mode volume 
[2"9"] . As <?jv — > Jj, losses due to energy relaxation and 
coupling to subradiant states will dominate the dynamics 
and superradiant effects will become more difficult to ob- 
serve. However as the intensity scales as, I^ ax ex N 2 , the 
resulting dynamics will be the determined from the com- 
petition between the superradiant intensity enhancement 
and the aforementioned losses. 

The proposal presented here is superior to existing 
superradiance experiments in three key aspects: initial 
state preparation, measurement and losses. In circuit- 
QED arbitrary initial states can be prepared at high 
fidelity using a sequence of quantum logic gates [17], 
whereas, in bulk solids and atomic gases a pulse is fired 
into the medium, randomly exciting sections of the en- 
semble O H014T2] . The lack of control over the ini- 
tial state in the latter schemes severely limits the pos- 
sible initial states that can be prepared and can lead 
to light being re-absorbed by the atomic system result- 
ing in delays and distortions of the exiting superradi- 
ant pulse, making analysis extremely complicated. Re- 
cent advances in circuit-QED have also lead to high fi- 
delity, single-shot joint multi-qubit measurements |17j . 
which can be used to completely characterise the state of 
the ensemble during emission. This can also be used to 
monitor the build up of field mediated correlations be- 
tween individual qubits. Such quantum correlations are 
a hallmark of the superradiant process [2 [3] . In atomic 
demonstrations of superradiance, characterisation of the 
quantum states of the atomic system using techniques 
such as quantum state tomography has not been possi- 
ble. Instead one monitored the emission from the ensem- 
ble to characterize superradiance. However in the su- 
perconducting case, full quantum state tomography can 



be performed on the superconducting qubits and detailed 
superradiance mediated correlations can be fully mapped 
out. Lastly, circuit-QED systems have relatively uniform 
qubit-field coupling rates and remarkably low rates of de- 
coherence and energy relaxation. In other implementa- 
tions, atomic collisions combined with spatially varying 
qubit-field coupling rates lead to losses which can only 
be overcome by using large ensembles. By using the low 
loss circuit-QED system presented here, it is no longer 
necessary to use such large ensembles to observe super- 
radiance. Small sample superradiance may have applica- 
tions in quantum information, including decoherence free 
subspace state encoding using subradiant states [5] and 
superradiant state readout [9]. 




Discussion of the two level approximation to 
the transmon 
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FIG. 3: Intensity of radiation from the resonator l3(Rit) ob- 
tained from numerical solution of the complete master equa- 
tion with the two level and three level approximations for each 
transmon. Parameters used are g^^/2-K — (83.7,85.7,85.1) 
MHz [19] and (u r , a r , k, 7J, 7J)/2tt = (6000, 660, 2000, 0.19, 0) 
MHz. The coupling rate to the third transmon level, Gz,j, is 
approximated using G3J = V2g3j |37j . 



In this paper we make the two level approximation 
to the transmon qubit, whereas infact the transmon is 
an anharmonic oscillator [30]. With several excitations 
in the system it may be possible to induce transitions to 
higher energy levels than the qubit levels, slowing the su- 
perradiant decay. Here, we analyse the effect that a third 
level for each transmon has on the dynamics. Firstly, we 
note that the two level approximation is a very common 
approximation in the literature [3TJH3TJ and has been ver- 
ified in several experiments [13 G23 EH]- Secondly, the 
transmon's anharmonicity can be varied considerably by 
altering the ratio of the Josephson energy to the charging 
energy (Ej/Ec)- From [3D] it can be seen that the anhar- 
monicity, a r , can be as large as 11% of the qubit level's 
transition frequency for Ej/Ec ~ 17.5 (i.e a r /27r ~ 660 
MHz at a qubit transition frequency of uj r /2n = 6GHz). 
This can be achieved with a modest increase in the trans- 
mons susceptibility to charge noise. At this level of an- 
harmonicity, the third level is far detuned from the tran- 
sition frequency and has very little effect on the system 
dynamics. 

In the bad cavity limit the photons escape immediately 
after the transmon decays, therefore there is no possibil- 
ity for the decay of one transmon to excite higher lev- 
els of another transmon. In our simulations, k is much 
larger than g^ to ensure this limit is valid. However, 
as k is finite, there is a small probability of the decay 
of one transmon exciting another transmon to a higher 
level. To see the effect, third level transmons in our sim- 
ulations. We hnd that the population of the third level 
is negligible. However, we do find the third level induces 
small oscillations in the intensity from the cavity (Fig. [3]) . 
This is a well known effect [35] [36] whereby a far detuned 
three level system acts like a two level system with minor 
oscillations at a frequency proportional to the detuning 
of the third level. 



III. PHASE MULTISTABILITY 

We now study another collective effect that can be ob- 
served in the same system - phase multistability. Phase 
multistability is the phenomenon where a coherently 
driven multi-qubit cavity system displaces the cavity field 
in phase depending on the collective state of the qubits. 
To observe phase multistability we return to the g > k 
regime in Eq. (nl) and resonantly drive the system. The 
system is now described by the following master equa- 
tion, 

N 

P = -i(529N,j(vJa ! + aaf),p] + [£(a f -a),p] (9) 

N s 

+ 1-D[a]p + J2-fV[cr-]p 

j=i 

where, we assume the drive field and the qubits are reso- 
nant with the TLR (A r , A gj = 0). Furthermore, we note 
that dephasing is negligible for transmon qubits, 7^ ~ 0. 

Single qubit phase bistability has proven very diffi- 
cult to experimentally achieve as the qubit-cavity sys- 
tem needs to fulfill the requirements of the strong cou- 
pling regime, g > k, 7 s . To date only a single experiment 
has detected signatures of phase bistability [2J, however, 
the modest ratio of gj k and various other experimental 
difficulties meant that the results were only marginally 
conclusive. Phase multistability for N > 1 has yet to be 
observed. In this section we show that the circuit-QED 
system presented here, with the addition of coherent driv- 
ing can display phase multistability and should present 
very clear experimental signatures. 



A. Phase Bistability 

Phase bistability can be described in the 7 s — >• 
limit as follows: at large intra-cavity photon numbers, 
n, the difference in energy between successive Jaynes- 
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Cummings manifolds \n, ±) f-> |n+l, ±) is E n+1 — E u 
uj r ± g/2^Jn, where w r is the frequency of the resonator 
and \n,±) — l/V2(\n + l,g) ± \n, e)). Strong driving 
of the TLR at tu r will populate large n Jaynes Cum- 
mings eigenstates and cause transitions along two sepa- 
rate pathways, \n, ±) -n- \n + 1,±). This driving is de- 
tuned from the transition frequency of each path by the 
manifold dependent detuning, ±g / (2yfn) . As the field is 
displaced in phase depending on the detuning, the field 
acquires a phase depending on the state of the qubit. The 
steady state density matrix of the system (p ss ) can be 
approximated by a mixture of two uncoupled, damped 
harmonic oscillators. These oscillators are driven off- 
resonantly and displaced in phase due to the different 
detunings of the effective driving fields experienced by 
each oscillator [T3] 14 . As t —> 00 the coherent state 
amplitude of each oscillator |aj,) is af s — f{2£f±ig)/n 1 
where / = \Jl — (g/2£) 2 . The two coherent states \af s ) 
can be detected by homodyne measurement of the res- 
onator field [ggPO] . 

Now relaxing the assumption, 7 s — > 0, we consider 
the case where, / — > 1, and in the steady state the 
system is in |ct+,+). The system will remain in this 
state until energy relaxation causes the qubit transition 

|+) — ► l/\/2(|+) + |— ))• After a short transient regime, 
the system will either reach the steady state, \af s , +) or 
|a~,— ) with equal probability. Due to this effect, in a 
single quantum trajectory the imaginary part of the cav- 
ity amplitude randomly switches between ±<?/«; at the 
rate Y/^M- 

To observe phase bistability experimentally, it is neces- 
sary to resolve the two coherent states \af s ). For strong 
driving this requires the ratio g/n to be as large as possi- 
ble. Also, as energy relaxation couples the two transition 
pathways, it can be shown |15j that phase bistability only 
occurs when 7 s < 2k. As circuit-QED systems can ful- 
fil each of these requirements, and does not suffer from 
problems associated with moving atoms in cavity QED 
[4T] . circuit-QED is a preferable system for the observa- 
tion of phase bistability. 

For large driving amplitude the steady state Q- 
distribution can be found after appropriate transforma- 
tions of the density matrix [13 , 



Q(x + iy) = 
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2^ir/3(£,£).U(l-*»)Ci-fc> 



dz, (10) 



where, Q{a) = {a\p as \(x) /V, a = x + iy and j3(a, b) is the 
beta function. A cross section of this function is shown 
in Fig. E] for realistic circuit-QED parameters. It is clear 
that \ocjg) can be resolved using homodyne measurement 
of the field [5S1 HT| . In contrast to |U| , we see that a 




FIG. 4: The Q-function in Eq. ( 10 1 at 35 = 2£/k with pa- 
rameters (g,K,j s )/2ir = (85,4, 0.19) MHz. Inset: The loca- 
tion of the two peaks in the steady state Q-function for a 
range of driving strengths £/2ir — 45 — 117.6 MHz, increas- 
ing from bottom to top. The following parameters were used: 
(3,k,7 s )/2tt = (85,28.3,0.19) MHz. 



circuit-QED demonstration of this effect will be dramatic 
and unambiguous demonstration of phase bistability. 

The inset in Fig. p] compares the peak locations af s 
to those obtained from numerical solution of the steady 
state density matrix for a range of driving amplitudes. At 
small driving amplitudes the system is too anharmonic 
to approximate as two harmonic oscillators. However, 
for larger amplitudes the two solutions coincide. Phase 
bistability has been the basis for several proposals in- 
cluding ultralow energy optical switching [41] . quantum 
feedback and qubit measurement [4^] . 



Similar to section II A we review the validity of the two 
level approximation to the transmon. The higher levels 
will lead to more transition pathways and therefore more 
peaks in the steady state Q-distribution. However, if the 
transmon is strongly anharmonic, these transitions will 
be suppressed and the additional peaks will be very mi- 
nor. As strong driving can populate higher lying states, 
to ensure the validity of the two level approximation the 
drive frequency needs to be significantly less than the an- 
harmonicity, a r ^$> £ > g/2 > k > 7 s /2. Furthermore, 
the populations of the higher levels can be reduced by 
enhancing the rate of energy relaxation 7 s . 



B. Phase Multistability 

When there are N qubits in the resonator an analo- 
gous phenomenon occurs: phase multistability [15] . Sim- 
ilarly to the single qubit case, phase multistability results 
from the energy structure of Tavis-Cummings manifolds 
at large excitation. For a given I, there are, 21 + 1, tran- 
sitions between the Dicke states \n)\l,m) ■(-> \n+ l)\l,m) 
PP. The difference in energy between successive Tavis- 




FIG. 5: Steady state Q-function for the field of a TLR with 
three qubits obtained from numerical solution of equations 
M with He. Center peaks correspond to a sa ' , whilst 

the other peaks correspond to a ss ' (see Eq. (ilk). Pa- 

rameters used are (£, k,7|)/2tt = (169.6,42.4,0.19) MHz and 
g:i,j/2-ir = (83.7,85.7,85.1) MHz [T9]. Note that the distance 
between the peaks can be increased by choosing smaller k. 



Cummings manifolds at large intra-cavity photon num- 



ber is, E. 



(l,rti) j-,(l,m) 
-C/71 



before, 



ui r +mgN J '■y/n [TS]- Proceeding as 
assuming 7" — > and gpjj sa g N , strong resonant 
driving of the TLR leads to transitions on 2/ + 1 sepa- 
rate ladders. As the field is displaced in phase depending 
on the detuning, the field acquires a phase depending on 
the state of the qubit. The steady state density matrix 
of the system can be approximated by a mixture of 21 + 1 
damped, driven uncoupled harmonic oscillators. Each 
harmonic oscillator has the coherent state amplitude, 



M,m) _ 



2/m(^/m +img N )/K, 



(11) 



Due to the small ratio of g~N / 'k and the inability to 
couple several qubits identically to a common field mode, 
phase multistability has not been experimentally demon- 
strated. However, using existing circuit QED param- 



where, f m = y/l - (mg N /£) 2 



eters Fig. IIIB demonstrates that phase multistability 
can be observed. The four peaks in the steady state 
Q-distribution in Fig. |IHB| are the four field displace- 
ments in phase space corresponding to the collective 
qubit states with I = -3/2, -1/2, 1/2, 3/2. The peak po- 
sitions a S s m coincide with the numerical solution to high 
precision. The centre two peaks are larger because the 
m = 1/2 transitions are driven closer to resonance than 
the m = 3/2 transitions. Phase multistability may have 
similar applications as phase bistability, i.e in quantum 
feedback |42j and ultralow energy microwave switching 
[HI H3J 03]. Moreover, it may be possible to use phase 
multistability to perform a collective measurement of the 
system to determine the collective spin, m |42j . 



IV. CONCLUSION 

In summary, we have demonstrated that two previ- 
ously unobservable collective quantum optical phenom- 
ena, small sample superradiance and phase multistabil- 
ity, can be observed in a circuit-QED based system us- 
ing existing experimental parameters. The proposal pre- 
sented here should provide a definitive and unambiguous 
demonstration of superradiance, with the added benefit 
of full quantum state tomography on the qubits. Owing 
to its the experimental tractability we expect that this 
approach will be indispensable for a range experiments 
on collective phenomena. 
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